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We argue that a normal Fermi liquid and a singlet, spin gapped Mott insulator cannot be con- 
tinuously connected, and that some intermediate phase must intrude between them. By explicitly 
working out a case study where the singlet insulator is stabilized by orbital degeneracy and an 
inverted Hund's rule coupling, mimicking a Jahn- Teller effect, we find that the intermediate phase 
is a superconductor. 
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Understanding metal-insulator transitions driven by 
electron correlation, the so-called Mott transitions (MT), 
is a long standing problem. Renewed attention was re- 
cently aroused by the development of the so-called Dy- 
namical Mean Field Theory (DMFT) @], a quantum ana- 
log of classical mean field theories, that treats exactly 
local temporal fluctuations while freezing spatial corre- 
lations. The DMFT has provided us with a complete 
characterization of the MT for the single band Hubbard 
model (SBHM) within the paramagnetic sector, where 
magnetic instabilities are not allowed. The insulator- 
metal MT does not occur by gradual closing of the Mott 
gap between the Hubbard bands in the single-particle 
spectrum, but it is associated with the appearance (at a 
given value U — U C 2) of a narrow "Kondo" quasiparti- 
cle peak or resonance at the Fermi level ||,[| . The peak 
height is pinned to the non- interacting value po(no), but 
its width is finite in the metal and vanishes continuously 
as the MT is approached. For U < U C 2 the metallic 
state is stable with respect to the metastable insulating 
solution which exists for U > U c i [@J§- Roughly in the 
interval U c i < U < U C 2, the spectral function presents 
both the broad Hubbard bands and the quasiparticle res- 
onance, thus combining to some extent the properties of a 
narrow-band metal with those of the Mott insulator. For 
smaller U, the lower and upper Hubbard bands merge 
together swamping the resonance away, and the corre- 
lated metal continuously turns into an ordinary metal. 
As soon as the constraint to the paramagnetic subspace 
is released, antiferromagnetism (AFM) appears, and ev- 
erything changes. Even in partially frustrated models, 
where the effects of magnetism are attenuated and AFM 
usually sets in only above some Um < U C 2, magnetism 
still "contaminates" the MT. The characteristic energy 
scale of the magnetic fluctuations, in proximity of a con- 
tinuous MT, necessarily exceeds the width of the Kondo- 
resonance, thus affecting the properties of the ensuing 



metallic phase, e.g., turning it into a spin-density wave. 

A new and different situation is to be expected if one 
could realize a singlet, spin-gapped Mott insulator, such 
as, e.g., a spin- liquid insulator. The entropy of such a 
state is zero, as opposed to the extensive entropy of the 
Mott phase in the SBHM Q , and no symmetry, magnetic 
or other, is broken. In this Letter we discuss the conse- 
quences of the zero-entropy property on the MT and on 
the correlated metal just below the transition point. To 
this end, we consider a (threefold) orbitally degenerate 
model arising in the physics of the alkali-doped fullerides 
||. For densities (n) =2 or 4, the ground state (GS) 
of this model system is nonmagnetic both for small U 
(metallic), and for large U (singlet nondegenerate Mott 
insulator). We show that a third, intermediate phase ac- 
tually intrudes between the Fermi liquid (FL) metal and 
the singlet Mott insulator. Remarkably, the intermediate 
phase is a superconductor. 

The Hamiltonian reads 
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where rn a = d\ aa d iaa is the electron number on each 
orbital (a = 1,2,3) at site i, rii = nu + na + 1143, and 
Ai a b = J2<r(di aa d ibcr + H.c). Besides the overall on-site 
repulsion U, the model includes a multiplet-exchangc- 
splitting term J. In presence of a Jahn- Teller coupling 
to some high frequency on-site vibration, as shown in Ref. 
W, J may effectively change sign from positive (Hund's 
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rule) to negative. When (n) corresponds to 2 or 4 elec- 
trons per site (the problem is electron-hole symmetric 
around (n) — 3), the isolated ion has a non degenerate 
GS which is simultaneously a spin and an orbital singlet. 
With a negative effective J, it is quite natural that a MT 
from the paramagnetic metal towards a non degenerate 
singlet insulator should take place increasing U . 

We carry out a DMFT study of the dynamics of model 
(|l|) in proximity of the MT for (n) = 2, 4 on a Bcthc 
lattice with hopping diagonal in the orbital index, t"? = 
tS a p. The bandwidth is W = 4i. The DMFT maps (Q) 
onto a threefold-degenerate Anderson impurity model 



(Vylc ana + H.c.) 
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where Hint is the interaction (|^) for the impurity opera- 
tors d OCT 's. The self-consistency condition is 
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G a {iu>) being the impurity Green's function for orbital a 
\ 1| . We solve the impurity model by exact diagonalization 
[3|. The sums in Eqs.(^) and (Q) are truncated to a finite 
and small number N s — 1 of conduction "bath" orbitals 
|jj . A surprisingly small value of N s is indeed enough to 
capture the qualitative features of the MT in the SBHM, 
and N s =5, the value used here, gives good quantita- 
tive results Although the Hamiltonian is 0(3) invari- 
ant under rotations in orbital space, we cannot exclude 
a spontaneous breaking of this symmetry. We explicitly 
studied solutions with orbital symmetry breaking (orbital 
ordering) , but we found that the orbitally symmetric so- 
lution, G a (ioj n ) = G(iu> n ), has always lower energy. In 
the following, we only consider such symmetric solutions. 

In order to characterize the metallic phase close to 
the MT, we compute the single-particle spectral func- 
tion p(u>) = — 1/tt TmG(u>), and the quasiparticle weight 
Z = m/m* at the Fermi level, which is finite in the metal- 
lic phase but zero in the insulator, and thus determines 
the position of the MT. Within our numerical accuracy 
the vanishing of Z appears to be continuous (see Fig. ^ 
(a)), as in the SBHM ||. In Fig. [l] we show the evolution 
of the spectral function across the transition, for different 
values of U/W but at fixed ratio J/W = -0.02. 

The metallic phase close to the transition (U < U C 2) 
shows a coexistence between high-energy insulating fea- 
tures (Hubbard bands), and the low-energy metal fea- 
ture, the Kondo resonance, whose width vanishes at the 
transition. In addition, the Hubbard bands display a 
well defined multiplet structure, absent in the SBHM, 
and typical of the orbitally degenerate site. Surprisingly, 



the linewidth of these atomic-like excitations is not set 
by the bandwidth as in the SBHM §, but by a much 
smaller energy scale, as if a kind of motional line narrow- 
ing were effectively at work. The coexistence of atomic 
and metallic features is therefore much more striking for 
our orbitally degenerate model than for the SBHM. 




FIG. 1. U dependence of the single-particle spectral func- 
tion for (n) = 4 and J/W = —0.02 in the normal state (with 
gauge symmetry enforced). 

As in the SBHM, the single-particle density of states 
(DOS) at the chemical potential p must coincide with 
the bare DOS po (/•*)• Moreover, as a consequence of Lut- 
tinger's theorem, the chemical potential must coincide 
with the bare one, so that p(p) = po{po) for every U, and 
the MT occurs by narrowing of the resonance peak. The 
peak width E F ~ ZW/2, may be viewed as an effective 
Fermi energy for the interacting system, or a Kondo tem- 
perature Tk, below which the system may be described 
as a FL. As a consequence, the entropy S(T)/N ~ T/Tk 
for T < Tk, so that S{Tk) is of order one up to C/ C 2, 
where Tk vanishes. By continuity, one can expect this 
entropy to be released at the MT. This is indeed what 
happens in the SBHM, where the entropy of the param- 
agnetic insulator is S/N — log 2. On the other hand, this 
is incompatible with our non-degenerate insulator, which 
has zero entropy. Therefore, it is hard to understand how 
a FL solution could continuously connect to the Mott in- 
sulator, despite the numerical evidence strongly pointing 
towards a continuous transition. This suggests that the 
metallic solution may be metastable close to the MT, and 
that some new phase may appear. 

The above heuristic arguments, which apply in any di- 
mension for a Brinkman-Rice type transition, where the 
m* diverges, can be generalized to other cases where the 
MT occurs towards a zero-entropy insulator. The main 
inadequacy of the metallic solution apparently lies in the 
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constraint imposed by the Luttinger theorem. Without 
it, the chemical potential could gradually move towards 
the edge of the DOS, thus allowing to smoothly connect 
the metal with a zero-entropy insulator, as for a metal 
to a band insulator transition. This scenario would im- 
ply a break-down of the Luttinger theorem before the 
MT. If the Mott insulator would break some symmetry 
of the Hamiltonian, e.g., the spin SU(2) symmetry, the 
metallic phase close to the MT would likely break the 
same symmetry, allowing the transition to become of the 
metal-band insulator type. This is what happens in the 
Hubbard model on a two-dimensional triangular lattice 
[U, and also corresponds to what is generally observed 
experimentally |5). If, as in our system, the insulator 
does not break any symmetry, it is less clear what kind 
of instability to expect. The presence of orbital degrees 
of freedom could suggest a possible orbital ordering, but, 
as we mentioned before, solutions with broken orbital 
symmetry have always higher energies with respect to 
orbitally symmetric solutions, ruling out the possibility 
of orbital ordering. Moreover, being the insulator a spin 
singlet, there is no reason to expect a magnetic instability 
in the metal, as confirmed by the vanishing of the local 
spin susceptibility at the MT. 

The only symmetry left, which we have so far enforced 
in solving the self consistency equation, is the U(l) gauge 
symmetry. A deeper analysis of the metallic phase in the 
framework of Landau FL theory suggests that U(l) sym- 
metry might be spontaneously broken. In our orbitally 
degenerate model, the Landau functional must be invari- 
ant under spin SU(2) and orbital 0(3) symmetries. We 
define the local densities in momentum space 
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where a are Pauli matrices, and Li (i 
angular momentum operators in the L = 

d 2 ~ Y\,y and d$ 
troduce five additional densities = ^ d J aak Wi ia bd bak , 
i = 1, . . . , 5, where the matrices Wj's are v3(L^ — Ly), 
V3(L x L y + L y L x ), V3(L y L z + L z L y ), V~3(L Z L X + L X L Z ) 
and (L x + Ly — 2L 2 Z ), respectively, as well as their spin 

antisymmetric counterparts w^ k = J2 d lak W i,a,b^a(3d b/3k . 
In terms of these densities, the Landau functional reads 



E 



1 

2^ 

k,k' 



fkk> n kn k ' + fkk' &k ■ &k 



1 

+ 2 



^kk' h,k 



i=x,y,z k,k' 



1 



W l k> 



(5) 




FIG. 2. U dependence of various relevant quantities 
in the normal state (with gauge symmetry enforced) for 
J/U = —0.02. In panel (a) we show the quasiparticle weight 
Z, in panel (b) the local spin susceptibility x> m panel (c) 
the compressibility k, and in panel (d) the real part of the 
zero-frequency pairing susceptibility n(0). 

The local spin susceptibility is given by x/xo = 
m*/[m(l + Fq-)], where \o is the bare susceptibility, 
Fq — 6V pofQ-m^/m, /q being the I = Legendre trans- 
form of the Landau parameter. Since x actually vanishes 
at the transition, as shown in Fig. ^ (b), then the anti- 
symmetric Landau parameter F A is positive and diverges 
at the MT, probably as to*. This is surprising since, for 
small U > I J|, /<f = -Z7/3 - 10/9J is negative. All lo- 
cal orbital and spin orbital susceptibilities vanish at the 
transition, since the insulator has a total gap for spin 
and orbital excitations, so that the Landau parameters 
Hq, Hq, Gq and Gq all diverge from the positive side. 
In addition, similarly to the SBHM, the compressibility 
k/kq = m*/[rn(l + Ffi )], vanishes like to -1 at the tran- 
sition (Fig. H (c)). Hence Fq ~ also diverges. As a 
result, the (intra-orbital) quasiparticle scattering ampli- 
tude in the singlet-channel 
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becomes negative close to the MT, finally vanishing 
as (to*) -1 for U — > U C 2, while at weak coupling 
t ^> U 3> \J\, this quantity is positive (A s = 
U/2 + 2 J/3 > 0). This strongly suggests a pairing 
between quasiparticles. To check for this instability, 
we compute the dynamical pair susceptibility n(w) = 
(0|Aj (w — H)~ l A |0), using as pair wave- function the 
singlet GS of the atomic limit A = J2 a ^Oat^Oai - 

The static limit of the real part of n(u>), shown in 
Fig. I (d), chan ges sign in a narrow region just before 
the MT, turning from negative (stable metal) to positive 
(unstable). This is another clear indication of a super- 
conducting instability of the FL. To confirm completely 
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this hypothesis, we finally allow for a superconducting 
solution of the DMFT [Q. The lattice model is mapped 
onto an impurity coupled to a superconducting bath, or, 
equivalently, coupled to a normal bath through a normal 
and an anomalous term. 

Hf = E 6anC ana c an<r + (^ n( ^acr c ana + H.C.) 
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stabilize a solution with a finite value of any A n , then 
a superconducting solution is found. We find a stable 
superconducting GS in the interval 0.82 < U < 0.9 (for 
J/U = —0.02) in close proximity of the MT, as shown 
in Fig. |[ We find that, while the superconducting or- 
der parameter vanishes at the MT, the spin gap of the 
supercondutor, as computed from the local spin suscep- 
tibility, remains finite, and joins continuosly into that of 
the insulator, as shown in Fig. || Hence neither phases 
possesses spin entropy. A superconducting phase just 
before the MT is striking, as there is no sign of, or rea- 
son for, such an instability in the weak-coupling limit, 
where the intraorbital s-wave scattering amplitude is re- 
pulsive. The anti-adiabatic electron-phonon coupling im- 
plicit in the choice J < 0, being unretarded, should be 
overwhelmed by a large U ^> |J|. As the repulsion U 
increases, the bare electrons turn into quasiparticles. An 
effective overscreening of U takes then place, generating 
the attraction that causes the s-wave pairing. This is 
reminiscent of phenomena in other systems, such as the 
Hubbard chain with next-nearest neighbor hopping and 
the two leg Hubbard ladder . Generation of an attrac- 
tive interaction in the effective impurity model, which 
has an on-site mechanism of singlet formation competing 
with the Kondo screening, is similar to the two impurity 
Kondo model O). For reasons of space we must defer 
pursuit of these analogies to future work. The heuristic 
arguments presented earlier suggest that superconduc- 
tivity close to a singlet MT should be a phenomenon of 
wider generality than the simple model where we have un- 
covered it. Experimentally, in the (n) — 4 alkali fullcridc 
Rb4C6o, which can be transformed by a pressure of about 
10 Kbar from its narrow-gap, probably singlet Mott in- 
sulating state, to a metal, no superconductivity was 
reported so far down to 0.39 K [^2). Nonetheless, even 
if some features of our model (antiadiabatic Jahn- Teller 
effect, diagonal hopping) are quite unrealistic, we can ex- 
pect all the same a narrow superconducting pocket suf- 
ficiently close to the critical insulator-metal pressure, at 



least if the transition is second order, or not too strongly 
first order. The small superconducting pocket for (n) = 4 
might be connected to the much larger one for (n) =3, 
the well-known case of superconducting fullerides. 
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FIG. 3. Superconducting coupling A$c (solid line) and 
spin gap A 3p in (dashed line) as computed from the local spin 
susceptibility. The verical dashed lines mark the boundaries 
of the various phases. 
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